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Scilab numbering policy used in this document and the relation to the 
above book. 

Exa Example (Solved example) 

Eqn Equation (Particular equation of the above book) 

AP Appendix to Example (Scilab Code that is an Appednix to a particular 
Example of the above book) 

For example, Exa 3.51 means solved example 3.51 of this book. Sec 2.3 means 
a scilab code whose theory is explained in Section 2.3 of the book. 
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Chapter 4 

Electron Ballistics 


Scilab code Exa 4.1 Calculation of acceleration time taken and distance covered an 

clc ; clear ; 

//Example 4.1 

//Calculation of acceleration , time taken , distance 
covered and kinetic energy of an accelerating 
proton 

//given values 

m = l . 67*10"-27 ; //mass of proton in kg 

q = l . 602*10"-19 ;// charge of proton in Coulomb 

vl=0;// initial velocity in m/s 

v2=2.5*10"6; //final velocity in m/s 

E = 500; //electric field strength in V/m 

/ / calculation 

a=E*q/m; //acceleration 

d i s p (a , ’ acceleration of proton in (m/s"2) is:’); 

t = v2/a; // time 

disp(t, ’time (in s) taken by proton to reach the 
final velocity is: ’); 
x = a*t"2/2; // distance 

disp(x,’ distance (in m) covered by proton in this 
time is: ’) ; 
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KE=E*q*x; //kinetic energy 

disp (KE , ’ kinetic energy(in J) at the time is:’); 


Scilab code Exa 4.2 electrostatic deflection 

clc ; clear ; 

//Example 4.2 

/ / electrostatic deflection 
//given values 

Vl=2000;//in volts , potential difference through 
which electron beam is accelerated 
1=.04; //length of rectangular plates 
d=.015;// distance between plates 
V = 50;// potential difference between plates 
// calculations 

alpha = atan(l*V/(2*d*Vl))*(180/°/ 0 pi);//in degrees 
disp (alpha , ’ angle of deflection of electron beam is: 

v=5.93*10~5*sqrt(Vl); //horizontal velocity in m/s 
t=l/v;//in s 

disp(t,’ transit time through electric field is:’) 


Scilab code Exa 4.3 electron projected at an angle into a uniform electric field 

clc ; clear ; 

//Example 4.3 

//electron projected at an angle into a uniform 
electric field 
//given values 

vl=4.5*10~5;// initial speed in m/s 

alpha = 37*°/ 0 pi/180;// angle of projection in degrees 

E = 200; / /electric field intensity in N/C 

e = 1.6*10"-19; //in C 
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m=9.1*10~-31;//in kg 
a=e*E/m;// acceleration in m/s "2 

t =2* vl * sin ( alpha)/a ; //time in s 

disp(t, ’time taken by electron to return to its 
initial level is: ’) 

H=(vl~2* sin (alpha)* sin (alpha))/(2*a);// height in m 
disp(H, ’maximum height reached by electron is: ’) 

s=(vl~ 2 ) *(2*sin(alpha)*cos (alpha))/(2*a); // 

displacement in m 

disp(s,’ horizontal displacement (in m) when it reaches 
maximum height is:’) 


Scilab code Exa 4.4 motion of an electron in a uniform magnetic field 

clc ; clear ; 

//Example 4.4 

//motion of an electron in a uniform magnetic field 
//given values 

V = 200;// potent ial difference through which electron 
is accelerated in volts 
B = 0.01 ;// magnetic field in wb / m' 2 
e = 1.6*10"-19; //in C 
m=9.1*10~-31;//in kg 

v=sqrt(2*e*V/m);// electron velocity in m/s 
d i s p (v , ’electron velocity is:’) 
r=m*v/(e*B);//in m 

disp (r,’radius of path (in m)is:’) 


Scilab code Exa 4.5 motion of an electron in a uniform magnetic field acting at an 

clc ; clear ; 

//Example 4.5 
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//motion of an electron in a uniform magnetic field 
acting at an angle 
//given values 
v=3*10"7; //e1ectron speed 
B= . 23 ;// magnetic field in wb/m"2 

q=45* °/ 0 pi /180 ; //in degrees , angle in which electron 
enter field 
e = 1.6*10"-19; //in C 
m=9.1*10~-31;//in kg 
R=m*v* sin (q)/(e*B); //in m 
disp(R, ’radius of helical path is:’) 
p = 2*°/ 0 pi*m*v*cos(q)/(e*B) ; //in m 
disp(p,’ pitch of helical path (in m) is:’) 


Scilab code Exa 4.6 Magnetostatic deflection 

clc ; clear ; 

//Example 4.6 

//Magnetostatic deflection 
//given values 
D = . 03; / /deflection in m 
m=9.1*10~-31;//in kg 
e = 1.6*10"-19; //in C 

L=.15;// distance between CRT and anode in m 

l=L/2; 

V=2000; //in voltsin wb/ 

B = D* sqrt (2*m*V) /(L*l*sqrt (e)); //in wb/m"2 
d i s p (B , ’transverse magnetic field acting (in wb / m' 2) 
is : ’) 


Scilab code Exa 4.7 electric and magnetic fields in crossed configuration 
clc ; clear ; 
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//Example 4.7 

//electric and magnetic fields in crossed 
configuration 
//given values 

B = 2* 10 ~-3 ; //magnetic field in wb / rn' 2 
E = 3.4* 10 ~4 ; //e 1 e c t r i c field in V/m 
m=9.1*10~-31;//in kg 
e=l.6*10"-19; //in C 
v = E/B;//in m/s 

d i s p (v , ’ electron speed is:’) 

R=m*v/(e*B); //in m 

disp (R,’radius of circular path (in m) when electric 
field is switched off’) 
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Chapter 5 

Electron Optics 


Scilab code Exa 5.1 Electron refraction calculation of potential difference 

clc ; clear ; 

//Example 5.1 

//Electron refraction , calculation of potential 
difference 

//given values 

VI =250 ; //potentia 1 by which electrons are 
accelerated in Volts 

alphal =50* °/„pi /180 ; //in degree 
alpha2 =30* °/„pi /180 ; //in degree 
b=sin(alphal)/sin(alpha2) ; 

/ / calculation 
V2=(b~2)*V1; 
a=V2-VI; 

d i s p (a , ’ potential difference (in volts) is:’); 


Scilab code Exa 5.2 Cyclotron 
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clc ; clear ; 

//Example 5.2&5.3 

//Cyclotron, calculation of magnetic induction , 
maximum energy 

//given values 

f=12*(10“6);// oscillator frequency in Hertz 
r =. 53;// radius of the dee in metre 
q = l . 6*10"-19 ; //Deuteron charge in C 
m = 3.34*10" -27 ; //mass of deuteron in kg 
/ / calculation 
B = 2*°/ 0 pi*f*m/q; // 

disp (B, ’ magnetic induction (in Tesla) is:’); 
E=B"2*q"2*r"2/(2*m); 

disp(E, ’maximum energy to which deuterons can be 
accelerated (in J) is’) 

E1 = E*6.24*10"18/10"6; //conversion of energy into MeV 
disp(El, ’maximum energy to which deuterons can be 
accelerated (in MeV) is’); 


Scilab code Exa 5.4 calculation of linear separation of lines formed on photograph 

clc ; clear ; 

//Example 5.4 

//Mass spectrograph , calculation of linear 

separation of lines formed on photographic plates 

//given values 

E=8*10"4; //c1ectric field in V/m 
B= . 55// magnetic induction in Wb/m*2 
q=l.6*10"-19; //charge of ions 

ml=20*1.67*10~-27; //atomic mass of an isotope of 
neon 

m2 = 22* 1.67* 10~-27 ; //atomic mass of other isotope of 
neon 
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11 //calculation 

12 x=2*E*(m2-ml)/(q*B~2); // 

13 disp(x, ’separation of lines (in metre) is:’) 
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Chapter 6 

Properties of Light 


Scilab code Exa 6.1 Optical path calculation 

clc ; clear ; 

//Example 6.1 

//Optical path calculation 
//given values 

n = l. 33; //refractive index of medium 
x=.75; //geometrical path in micrometre 
/ / calculation 
y=x*n; // 

di sp (y , ’ o p t i c al path (in micrometre) is:’) 


Scilab code Exa 6.2 Coherence length calculation 

clc ; clear ; 

//Example 6.2 

//Coherence length calculation 
//given values 
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6 1 = 1*10"-14;// line width in metre 

7 x = 10.6*10“-6; //IR emission wavelength in metre 

8 //calculation 

9 y=x"2/1; // 

10 disp (y coherence length (in metre) is:’) 
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Chapter 7 

Interference and Diffraction 


Scilab code Exa 7.1 plane parallel thin film 

clc ; clear ; 

//Example 7.1 

//plane parallel thin film 
//given values 

x=5890*10~-10; //wavelength of light in metre 
n = 1.5; / /refractive index 

r = 60* %pi /180 ;// angle of refraction in degree 
/ / calculation 

t = x/ (2*n*cos (r)) ; 

disp(t*10~6,’ thickness of plate (in micrometre) is:’ 

) ; 


Scilab code Exa 7.2 wedge shaped thin film 

clc ; clear ; 

//Example 7.2 

//wedge shaped thin film 
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//given values 

x=5893*10~-10 ;// wavelength of light in metre 
n = 1.5; / /refractive index 
y = . 1 * 10 ~ - 3 ; / / fringe spacing 
/ / calculation 

z=x/(2*n*y);// angle of wedge 

alpha = z* 180/°/ 0 pi ;// conversion of radian into degree 
disp(alpha, ’angle of wedge (in degree) is:’); 


Scilab code Exa 7.3 Newtons ring experiment 

clc ; clear ; 

//Example 7.3 

//Newton’s ring experiment— calculation of 
refractive index 

//given values 

D1 = 1.5 ; //diametre (in cm) of tenth dark ring in air 
D2 = 1.27 ; //diametre (in cm) of tenth dark ring in 
liquid 


/ / calculation 

n = D1~2/D2 ~2; 

di sp (n , ’ r e f r a c t i v e index of liquid is’); 


Scilab code Exa 7.4 nonreflecting film 

clc ; clear ; 

//Example 7.4 
/ / nonreflecting film 
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5 //given values 

6 1 = 5500*10"-10;// wavelength of light 

7 nl = 1.33 ; //r e fr ac t i ve index of water 

8 n2=1.52; //refractive index of glass window pane 

9 x = sqrt(nl);//to check if it is nonreflecting 
10 

11 //calculation 

12 t = l/(4*nl);// thickness of water film required 

13 disp(t*10~6, ’minimum thickness of film (in metre) 

’); 


i s 
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Chapter 8 
Polarization 


Scilab code Exa 8.2 Polarizer 

clc ; clear ; 

// Example 8.2 

//Polarizer , calculation of angle 
//given values 

10 = 1; //intensity of polarised light 

11 = Io/2; // intensity of beam polarised by first by 

first polariser 

12 = Io/3 ; //i n t e nsi t y of light polarised by second 

polariser 


/ / calculation 

a=acos (sqrt (12/11)); 

alpha = a* 180/°/ 0 pi ;// conversion of angle into degree 
disp (alpha , ’ angle between characteristic directions 
(in degree) is’); 


Scilab code Exa 8.3 calculation of birefringence 
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clc ; clear ; 

//Example 8.3 

//calculation of birefringence 
//given values 

1 = 6* 10~-7 ;// wavelength of light in metre 
d=3*10~-5; //thickness of crystal 


/ / calculation 

x = 1/(4*d) ; 

disp(x,’the birefringance of the crystal is’); 
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Chapter 11 

Architectural Acoustics 


Scilab code Exa 11.1 calculation of total absorption and average absorption coeffl 

clc ; clear ; 

//Example 11.1 

//calculation of total absorption and average 
absorption coefficient 

//given values 

V = 20* 15*5 ;// volume of hall in nT3 

t=3.5;// reverberation time of empty hall in sec 


/ / calculation 

al = .161*V/t; //total absorption of empty hall 
k=al/(2*(20*15+15*5+20*5)); 

disp(k,’the average absorption coefficient is’); 


Scilab code Exa 11.2 calculation of average absorption coefficient 
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clc ; clear ; 

//Example 11.2 

/ / calculation of average absorption coefficient 
//given values 

V = 10*8*6 ;// volume of hall in m'3 
t = l. 5; //reverberation time of empty hall in sec 
A = 20;//area of curtain cloth in m' 2 
tl = l;//new reverberation time in sec 

/ / calculation 

al = .161*V/t; //total absorption of empty hall 
a2=.161*V/tl; //total absorption after a curtain 
cloth is suspended 

k=(a2-al)/(2*20); 

disp(k,’the average absorption coefficient is’); 


Scilab code Exa 11.3 calculation of average absorption coefficient and area 

clc ; clear ; 

//Example 11.3 

//calculation of average absorption coefficient and 
area 

//given values 

V = 20*15*10 ; //volume of hall in m'3 

t=3.5;// reverberation time of empty hall in sec 

11 = 2.5 ; //reduced reverberation time 

k2=. 5; //absorption coefficient of curtain cloth 

II calculation 

al = .161*V/t; //total absorption of empty hall 

kl = al/(2*(20*15+15*10 + 20*10) ) ; 
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14 disp(kl,’the average absorption coefficient is’); 

15 a2=.161*V/tl; //total absorption when wall is covered 

with curtain 

16 a=t 1 * ( a2 -al) / (11 *k2) ; 

17 disp(a, ’area of wall to be covered with curtain (in m 

'2)is : ’) 
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Chapter 12 
Ultrasonics 


Scilab code Exa 12.1 calculation of natural frequency 

clc ; clear ; 

//Example 12.1 

//calculation of natural frequency , magnetostriction 
//given values 

1=40*10“-3;// length of pure iron rod 
d = 7.25*10~3; //density of iron in kg/m“3 
Y = 115* 10“9 ; //Young ’ s modulus in N/m' 2 

/ / calculation 

f = (l*sqrt (Y/d))/(2*1) ; 

disp(f*10“-3, ’the natural frequency (in kHz) is’); 


Scilab code Exa 12.2 calculation of natural frequency 


clc ; clear ; 

/ / Example 12.2 
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/ / calculation of natural frequency 
//given values 

t=5.5*10~-3; //thickness in m 
d=2.65*10~3; //density in kg/m~3 
Y = 8*10" 10 ; //Young ’ s modulus in N/rn' 2 


/ / calculation 

f = (sqrt(Y/d))/(2*t);//frequency in hertz 
disp(f*10~-3, ’the natural frequency (in kHz) is’); 


Scilab code Exa 12.3 calculation of depth and wavelength 

clc ; clear ; 

/ / Example 12.3 

//calculation of depth and wavelength 
//given values 

f = .07*10~6; / / frequency in Hz 

t=. 65; //time taken for pulse to return 

v = 1700;// velocity of sound in sea water in m/s 

/ / calculation 
d=v*t/2; // 

disp(d,’the depth of sea (in m) is’); 
l = v/ f;//wavelenght of pulse in m 
disp(l*10~2, ’wavelength of pulse (in cm) is’); 
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Chapter 13 
Atomic Physics 


Scilab code Exa 13.1 calculation of rate of flow of photons 

clc ; clear ; 

//Example 13.1 

//calculation of rate of flow of photons 
//given values 

1 =5893* 10~ -10 ;// wavelength of light in m 
P = 40;//power of sodium lamp in W 
d = 10;// distance from the source in m 
s=4*°/„pi *d~2 ; // surface area of radius in m' 2 
c=3*10~8; //velocity of light in m/s 
h = 6.626*10~-34 ; //Planck ’ s constant in Js 
/ / calculation 
E = P * 1; // 

di sp (E , ’ t o t a 1 energy emitted per second (in Joule) is’ 

) ; 

n = E*1/(c*h) ; //tota 1 no of photons 

R=n/s; 

disp(R, ’rate of flow of photons per unit area (in /m 
* 2) is’) 
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Scilab code Exa 13.2 calculation of threshold wavelength and stopping potential 


clc ; clear ; 

//Example 13.2 

//calculation of threshold wavelength and stopping 
potential 

//given values 

1=2000;// wavelength of light in armstrong 
e=1.6*10~-19; //charge of electron 
W = 4.2; //work function in eV 
c = 3*10~8; //velocity of light in m/s 
h = 6.626*10~-34 ; //Planck ’ s constant in Js 
/ / calculation 

x = 12400 / (W) ; //h*c = 12400 eV 

disp (x , ’ threshold wavelength (in Armstrong ) i s ’) ; 

Vs = (12400/1-W) ; // 

disp(Vs, ’stopping potential (in VOLTS) is ’) 


Scilab code Exa 13.3 calculation of momentum of Xray photon undergoing scattering 

clc ; clear ; 

//Example 13.3 

//calculation of momentum of X—ray photon undergoing 
scattering 

//given values 

alpha=60*%pi/180;// scattering angle in radian 
e=1.6*10~-19; //charge ofelectrone 
W = 12273 ;// work function in eV 
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c=3*10“8; / /velocity of light in m/s 
h=6.626*10“-34; //Planck’s constant in Js 
he = 12400; //in eV 

m = 9.1*10~-31// restmass of photon in kg 
/ / calculation 

x = hc/ (W) ;// wavelength of photon undergoing modofied 
scattering in arm strong 

y = x-(h/(m* c)) *(l-cos (alpha)); 
p = h/y* 10 “10; 

disp(p, ’momentum of photon (in kg—m/ s ) is’); 


Scilab code Exa 13.4 calculation of wavelength of scattered radiation and velocity 

clc ; clear ; 

//Example 13.4 

//calculation of wavelength of scattered radiation 
and velocity of recoiled electron 

//given values 

alpha=30*°/ 0 pi/180 ;// scattering angle in radian 
e=1.6*10“-19; //charge ofelectron 

x = l . 372*10“-10 ; //wavelength of incident radiation in 
m 

c=3*10“8; //velocity of light in m/s 
h = 6.626*10“-34 ; //Planck ’ s constant in Js 
m = 9 . 1*10“-31// rest mass of photon in kg 
he = 12400; //in eV 
/ / calculation 

y=((x+(h/(m* c)) *(l-cos (alpha))))* 10“10; 
disp (y , ’ wavelength of scattered radiation(in 
armstrong ) is ’) ; 

xl=x*10~10; //converting incident wavelength into 
arm strong 
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KE=hc*e*((1/xl)-(1/y)); //kinetic energy in Joule 
disp(KE,’ kinetic energy in joule is ’); 
v=sqrt (2*KE/m); 

disp(v,’ velocity of recoiled electron (in m/s"2)is 7 ) 


Scilab code Exa 13.5 calculation of wavelength of light emitted 

clc ; clear ; 

//Example 13.5 

//calculation of wavelength of light emitted 
//given values 

e=1.6*10~-19; //charge of electrone 
c=3*10~8;// velocity of light 
h=6.626*10"-34; //Planck’s constant in Js 
El=5.36; //energy of first state in eV 
E2=3.45; //energy of second state in eV 


/ /1) calculation 

l = h*c*10~10/((El-E2)*e) ; 

disp (1, ’ wavelength of scattered light (in Armstrong) 
is ’) ; 


Scilab code Exa 13.6 calculation of de Broglie wavelength 

clc ; clear ; 

//Example 13.6 

//calculation of de Broglie wavelength 
/ / 1)giv e n values 
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e = 1.6*10"-19; 

h = 6.626*10"-34 ; //Planck ’ s constant in Js 
V = 182; //potential difference in volts 
m = 9.1*10"-31 ; //mass of e in kg 


// 1) calculation 

l=h/ sqrt (2*e*m*V); 

d i s p (1 , ’ de Brogliewavelength (in m) i s ’ ) ; 


/ /2)given values 

ml = l;//mass of object in kg 

v = l; //velocity of object in m/s 

11 = h/(ml *v) ; 

disp (11 , ’ debrogie wavelength of object in m) is’); 


Scilab code Exa 13.7 calculation of uncertainty in position 

clc ; clear ; 

//Example 13.7 

//calculation of uncertainty in position 
//l) given values 

h = 6.626*10"-34 ; //Planck ’ s constant in Js 
vl=220; //velocity of e in m/s 
m = 9 . 1 * 10"-31 ; //mass of e in kg 
A=0.065/100; //accuracy 


II 1) calculation 

v2 = vl*A;// uncertainty in speed 
xl = h/ (2*°/ 0 pi*m*v2) ; / / 


34 



16 

17 

18 

19 

20 

21 

22 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 


disp (xl uncertainty in position of e (in m)is’); 


/ / 2) g i v e n values 

ml = 150/1000 ; //mass of object in kg 

x2 = h/ (2* °/ 0 pi *ml * v2) ; 

disp (x2 , ’ uncertainty in position of baseball (in m) 
is ’) ; 


Scilab code Exa 13.8 Energy states of an electron and grain of dust 

clc ; clear ; 

//Example 13.8 

//calculation of energy states of an electron and 
grain of dust and comparing 

//l) given values 

Ll = 10*10~-10;// width of potential well in which e is 
confined 

L2= . 1*10“-3; //width of potential well in which grain 
of dust is confined 

h=6.626*10"-34; //Planck’s constant in Js 
vl=10“6; //velocity of garin of dust in m/s 
ml =9.1 * 10" -31; //mass of e in kg 
m2 = 10"-9 ; //mass of grain in kg 

/ /1) calculation 

Ee1 = 1"2*h"2/(8*ml * LI~2) ; //first energy state of 
electron 

disp (Ee1 , ’ first energy state of e is ’); 

Ee2=2"2*h"2/(8*ml * LI"2) ; //second energy state of 
electron 

disp (Ee2second energy state of e is ’); 

Ee3 = 3"2*h"2/(8*ml * LI~2) ; //third energy state of 
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electron 

20 disp (Ee3third energy state of e is ’); 

21 disp( ’Energy levels of an electron in an infinite 

potential well are quantised and the energy 
difference between the successive levels is quite 
large . Electron cannot jump from one level to 
other on strength of thermal energy . Hence 
quantization of energy plays a significant role 
in case of electron’) 

22 

23 Egl = 1~2*h~2/(8*m2 *L2~2) ; //first energy state of 

grain of dust 

24 disp (Egl f i r s t energy state of grain of dust is ’); 

25 Eg2=2~2*h~2/(8*m2*L2~2); //second energy state of 

grain of dust 

26 disp(Eg2, ’second energy state of grain of dust is ’) 

> 

27 Eg3=3~2*h~2/(8*m2*L2~2);// third energy state of 

grain of dust 

28 disp(Eg3, ’third energy state of grain of dust is ’) 

> 

29 KE = ra2 * vl ~ 2/2 ; //k i n e t i c energy of grain of dust; 

30 disp(KE, ’kinetic energy of grain of dust is’); 

31 disp(’The energy levels of a grain of dust are so 

near to each other that they constitute a 
continuum . These energy levels are far smaller 
than the kinetic energy possessed by the grain of 
dust. It can move through all these energy levels 
without an external supply of energy . Thus 
quantization of energy levels is not at all 
significant in case of macroscopic bodies.’) 
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Chapter 14 
Lasers 


Scilab code Exa 14.1 calculation of intensity of laser beam 

clc ; clear ; 

//Example 14.1 

//calculation of intensity of laser beam 

//given values 
P = 10* 10~-3 ; //Power in Watt 
d = 1.3*10~-3;//diametre in m 
A = 0 /,pi *d ~ 2/4 ;// area in m' 2 


/ / calculation 

I=P/A; 

disp(I, ’intensity (in W/m'2) is’); 


Scilab code Exa 14.2 calculation of intensity of laser beam 

clc ; clear ; 

//Example 14.2 
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//calculation of intensity of laser beam 

//given values 
P=l*10~-3 ; //Power in Watt 
1=6328*10“-10 ;// wavelength in m 
A = l~2;//area in in' 2 


/ / calculation 

I=P/A; 

disp(I, ’intensity (in W/m'2) is’); 


Scilab code Exa 14.3 calculation of coherence length bandwidth and line width 

clc ; clear ; 

//Example 14.3 

//calculation of coherence length , bandwidth and line 
width 

//given values 

c = 3*10~8; //velocity of light in m/s 
t=.l*10“-9; //timedivision in s 
1=6238* 10“ -10 ;// wavelength in m 

/ / calculation 
x=c*t; 

disp(x, ’coherence length (in m) is’); 
d=l/t; 

disp(d, ’bandwidth (in Hz) is’); 

y=l~2*d/c; //line width in m 

disp(y*10“10,’line width (in arm strong ) i s ’) ; 


Scilab code Exa 14.4 calculation of frequency difference 
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clc ; clear ; 

//Example 14.4 

/ / calculation of frequency difference 
//given values 

c=3*10“8; //velocity of light in m/s 
1=.5; / /distance in m 

/ / calculation 
f=c/(2*l);//in hertz 

disp(f/10“6, ’frequency difference (in MHz) is ’) ; 


Scilab code Exa 14.5 calculation of frequency difference 

clc ; clear ; 

//Example 14.5 

//calculation of no of cavity modes 
//given values 

c=3*10“8; / /velocity of light in m/s 

n = 1.7 5; //refractive index 

1=2*10“-2;// length of ruby rod in m 

x = 6943* 10~-10 ;// wavelength in m 

y = 5.3*10“-10 ; //spread of wavelength in m 

/ / calculation 

d=c/n/1; 
f=c*y/x“2; 
m=f/d; 

disp (m, ’no of modes is’); 
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Chapter 15 

Atomic nucleus and nuclear 
energy 


Scilab code Exa 15.1 calculation of binding energy per nucleon 

clc ; clear ; 

/ / Example 15.1 

//calculation of binding energy per nucleon 
//given values 

Mp = 1.00814 ; //mass of proton in amu 
Mn = l . 008665 ; //mass of nucleon in amu 
M = 7.01822; //mass of Lithium nucleus in amu 
amu=931; //amu in MeV 

n = 7-3;//no of neutrons in lithium nucleus 
/ / calculation 

ET=(3*Mp+4*Mn-M)*amu; //total binding energy in MeV 
E = ET/7;//7 Is the mass number 

d i s p (E , ’Binding energy per nucleon in MeV is’); 
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Scilab code Exa 15.2 calculation of energy 

clc ; clear ; 

//Example 15.2 
//calculation of energy 

//given values 

Ml = 15.00001 ;// atomic mass of N15 in amu 
M2 = 15.0030 ;// atomic mass of 015 in amu 
M3 = 15.9949;// atomic mass of 016 in amu 
amu = 9 31.4; / / amu in MeV 
mp = l . 0072766 ;// restmass of proton 
mn = l . 0086654 ;// restmass of neutron 

/ / calculation 

Ql=(M3-mp-Ml)*amu; 

disp(Ql, ’energy required to remove one proton from 
016 is’); 

Q2=(M3-mn-M2)*amu; 

disp(Q2, ’energy required to remove one neutron from 
016 is’); 


Scilab code Exa 15.3 calculation of binding energy 

clc ; clear ; 

//Example 15.3 

//calculation of binding energy 
//given values 

Mp = l . 00758 ; //mass of proton in amu 
Mn = 1.00897 ; //mass of nucleon in amu 
M = 4.0028;// mass of Helium nucleus in amu 
amu = 9 31.4; / / amu in MeV 

II calculation 
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El=(2*Mp+2*Mn-M)*amu; //total binding energy 
disp(El, ’Binding energy in MeV is’); 

E2 = E1*10~6*1. 6*10“-19; 

disp (E2,’binding energy in Joule is’); 


Scilab code Exa 15.4 calculation of binding energy 

clc ; clear ; 

/ / Example 15.4 

//calculation of amount of unchanged material 

//given values 

T = 2; //half life in years 

k=.6931/T;// decay constant 

M = 4.0028;// mass of Helium nucleus in amu 

amu = 9 31.4 ; / / amu i n MeV 

No = l; // initial amount in g 

/ / calculation 

N = No*(%e''(-k*2*T)) ; 

disp (N, ’ amount of material remaining unchanged after 
four years (in gram) is’); 


Scilab code Exa 15.5 calculation of halflife 

clc ; clear ; 

/ / Example 15.5 

//calculation of amount of halflife 

//given values 
t=5;//time period in years 
amu = 9 31.4; / / amu in MeV 
No=5;// initial amount in g 
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N=5-(10.5*10~-3); //amount present after 5 years 


/ / calculation 

k=log(N/No)/t; //decay constant 

T=-.693/k; 

d i s p (T , ’ halflife in years is’); 


Scilab code Exa 15.6 calculation of activity 

clc ; clear ; 

//Example 15.6 
//calculation of activity 

//given values 

t=28;//half life in years 
m = 10~-3 ; //mass of sample 
M = 90;// atomic mass of strontium 
NA=6.02*10~26; //avogadro’ s number 


II calculation 

n=m*NA/M; //no of nuclei in 1 mg sample 
k=.693/(t*365*24*60*60); //decay constant 
A=k*n; 

disp (A , ’ act ivit y of sample ( in disintegrations per 
second ) is ’) ; 


Scilab code Exa 15.7 calculation of age of mineral 

clc ; clear ; 

/ / Example 15.7 

//calculation of age of mineral 
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//given values 

t=4.5*10~9 ; //half life in years 

Ml =238;// atomic mass of Uranium in g 

m=. 093;// mass of lead in 1 g of uranium in g 

NA=6.02*10~26; //avogadro’s number 

M2 =206;// atomic mass of lead in g 

/ / calculation 

n = NA/Ml;//no of nuclei in 1 g of uranium sample 
nl=m*NA/M2 ; //no of nuclei in m mass of lead 
c=nl/n; 

k=.693/t;// decay constant 

T=(l/k)*log(l + c) ; 

disp(T, ’age of mineral in years is ’); 


Scilab code Exa 15.8 calculation of age of wooden piece 

clc ; clear ; 

//Example 15.8 

//calculation of age of wooden piece 
//given values 

t=5730; //half life of C14 in years 
Ml=50;//mass of wooden piece in g 

Al=320; //activity of wooden piece (disintegration 
per minute per g) 

A2=12; //activity of living tree 

/ / calculation 

k=.693/t;// decay constant 

A=A1/M1;// activity after death 

T=(1/k) *log (A2/A); 

disp(T, ’age of mineral in years is ’) ; 
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Scilab code Exa 15.9 calculation of energy released 

clc ; clear ; 

//Example 15.9 

//calculation of energy released 
//given values 

Ml = 10.016125 ; //atomic mass of Boron in amu 
M2 = 13.007440 ;// atomic mass of C13 in amu 
M3 = 4.003874;// atomic mass of Helium in amu 
mp = l . 008146 ; //mass of proton in amu 
amu = 931; //amu in MeV 

/ / calculation 

Q=(Ml+M3-(M2+mp))*amu; //total binding energy in M 
d i s p (Q , ’ Binding energy per nucleon in MeV is’); 


Scilab code Exa 15.10 calculation of crossection 

clc ; clear ; 

/ / Example 15.10 
//calculation of crosssection 

//given values 

t=.01*10~-3; //thickness in m 

n = 10~13;//no of protons bombarding target per s 
NA=6.02*10~26; //avogadro’s number 
M = 7; //atomic mass of lithium in kg 
d = 500;// density of lithium in kg /rrh 3 
n0=10~8;//no of neutrons produced per s 
/ / calculation 
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nl = d*NA/M; //no of target nuclei per unit volume 
n2=nl*t;//no of target nuclei per area 
A=nO/(n*n2); 

d i s p (A , ’ crosssection (in rn' 2) for this reaction is’); 


Scilab code Exa 15.11 calculation of final energy 

clc ; clear ; 

/ / Example 15.11 
//calculation of final energy 

//given values 

B=.4;//max magnetic field in Wb/m' 2 

c=3*10~8; 

e = 1.6*10"-19; 

d = 1.52 ; // diametre in m 

r=d/2; 

/ / calculation 

E=B*e*r*c ; //E=pc , p=mv=Ber 

disp(E,’ final energy of e (in J) is’); 

El=(E/e)/10~6; 

disp(El, ’final energy of e (in MeV) is’); 


Scilab code Exa 15.12 calculation of amount of fuel 

clc ; clear ; 

//Example 15.12 

//calculation of amount of fuel 
//given values 

P = 100* 10~6 ; //power required by city 
M = 235;// atomic mass of Uranium in g 
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e = 20/100 ; / / conversion efficiency 
NA=6.02*10~26; //avogadros number 

E=200*10~6*1.6*10~-19; //energy released per fission 
t=8.64*10~4; //day in seconds 


/ / calculation 

El=P*t ; // energy requirement 

m = El*M/(NA*e*E); //no of nuclei N=NA*m/M, energy 

released by m kg is N*E, energy requirement=e*N*E 
di sp (m , ’ amount of fuel (in kg) required is’); 


Scilab code Exa 15.13 calculation of power output 

clc ; clear ; 

//Example 15.13 
//calculation of power output 

//given values 

M = 235;// atomic mass of Uranium in kg 
e = 5/100; / / reactor efficiency 

ra = 25/1000; //amount of uranium consumed per day in kg 
E=200*10~6*1.6*10~-19; //energy released per fission 
t=8.64*10~4; //day in seconds 
NA=6.02*10~26; //avogadros number 

/ / calculation 

n=NA*m/M; //no of nuclei in 25g 
El=n*E;// energy produced by n nuclei 
E2=E1*e;// energy converted to power 
P = E2/t ; //power output in Watt 
disp (P/10~6 , ’power output in MW is ’) ; 
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Scilab code Exa 15.14 calculation of power developed 


clc ; clear ; 

//Example 15.14 

//calculation of power developed 
//given values 

M = 235;// atomic mass of Uranium in kg 

ra = 20.4; //amount of uranium consumed per day in kg 

E=200*10~6*1.6*10~-19; //energy released per fission 

t =3600* 1000 ;// time of operation 

NA=6.02*10~26; //avogadros number 

/ / calculation 

n = NA*m/M; //no of nuclei in 20.4kg 
El=n*E; //energy produced by n nuclei 
P = El/t;//in Watt 

disp(P/10~6, ’power developed in MW is ’) ; 


Scilab code Exa 15.15 calculation of amount of dueterium consumed 

clc ; clear ; 

//Example 15.15 

//calculation of amount of dueterium consumed 
//given values 

Ml =2.01478 ;// atomic mass of Hydrogen in amu 

M2 =4.00388 ;// atomic mass of Helium in amu 

amu = 931; //amu in MeV 

e=30/100; / / efficiency 

P=50*10~6; //output power 

NA=6.026*10~26; //avogadro number 

t=8.64*10~4; //seconds in a day 

II calculation 
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15 Q = (2* Ml -M2) *amu ; //energy released in a D—D reaction 

i n MeV 

16 0 = Q*e*10~6/2;// actual output per dueterium atom in 

eV 

17 n = P/(0*l . 6*10“-19) ; //no of D atoms required 

18 m = n*Ml/NA; //equivalent mass of D required per s 

19 X=m*t ; 

20 

21 disp(X, ’Deuterium requirement per day in kg is ’); 
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Chapter 16 

Structure of Solids 


Scilab code Exa 16.1 calculation of density 

clc ; clear ; 

//Example 16.1 
//calculation of density 

//given values 

a=3.36*10~-10; //lattice constant in m 
M = 209 ; //atomicmass of polonium in kg 
N=6.02*10~26; //avogadro’s number 
z = l ; //no of atom 
/ / calculation 
d=z*M/(N*a~3) 

d i s p (d , ’ density (in kg/m "3) is’); 


Scilab code Exa 16.2 calculation of no of atoms 

clc ; clear ; 

//Example 16.2 
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//calculation of no of atoms 
//given values 

a = 4.3*10"-10 ; //edge of unit cell in m 
d = 963 ; // density in kg/m' 3 
M = 23; //atomicmass of sodium in kg 
N=6.02*10~26; //avogadro’s number 

/ / calculation 

z=d*N*a~3/M; 

disp (z, ’no of atoms is ’) ; 


Scilab code Exa 16.3 calculation of distance 

clc ; clear ; 

//Example 16.3 
//calculation of distance 

//given values 

z=4;//no of atoms in fee 

d = 2180 ;// density in kg/in' 3 

M = 23 + 35.3;// atomicmass of sodium chloride in kg 
N=6.02*10~26; //avogadro’s number 

/ / calculation 

al=z*M/(N*d); 
a=al~ ( 1/3) ; 

1 = a / 2 ; / / i n m 

di sp (1 * 10" 10 , ’distance between adjacent chlorine and 
sodium atoms in armstrong is ’); 
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Scilab code Exa 16.4 calculation of interatomic spacing 

clc ; clear ; 

//Example 16.4 

//calculation of interatomic spacing 
//given values 

alpha = 30*%pi/180 ;// Bragg angle in degree 

h = l; 
k = l; 

1 = 1 ; 

m = 1; //order of reflection 
x = l. 75* 10~-10 ;// wavelength in m 

/ / calculation 

d=m*x/ (2*sin(alpha)) ; 
a=d*sqrt(h~2+k~2+l~2) ; //in m 

disp(a*10~10, ’interatomic spacing in arm strong is ’) ; 
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Chapter IT 

The Band Theory of Solids 


Scilab code Exa 17.1 calculation of velocity of fraction of free electrone 

clc ; clear ; 

//Example 17.1 

//calculation of probability 
//given values 

E = . 01; / /energy difference in eV 

kT=. 026; //temperture equivalent at room temp in e 

/ / calculation 

P = l/(l+(°/„e~(E/kT))) ; 

disp (P , ’ interatomic spacing is’); 


Scilab code Exa 17.2 calculation of velocity of e 

clc ; clear ; 

//Example 17.2 

//calculation of velocity of e 
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//given values 

e=l.6*10~-19; //charge of e in C 
E=2.1*e;//fermi level in J 
m = 9.1*10~ -31 ; //mass of e in kg 

/ / calculation 

v=sqrt (2*E/m); 

disp(v,’ velocity of e(in m/s ) ’ ) ; 


Scilab code Exa 17.3 calculation of velocity of fraction of free electrones 

clc ; clear ; 

//Example 17.3 

//calculation of velocity of fraction of free 
electrons 

//given values 
E=5.5;//fermi level in eV 

kT=. 026; //temperture equivalent at room temp in e 

/ / calculation 

f=2*kT/E; 

disp (f , ’ fr act ion of free electrone\s upto width kT 
on either side of Ef is’); 
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Chapter 18 
Semiconductors 


Scilab code Exa 18.2 calculation of probability 

clc ; clear ; 

//Example 18.2 

//calculation of probability 

//given values 
T=300;//temp in K 

kT=. 026; //temperture equivalent at room temp in eV 
Eg=5.6; //forbidden gap in eV 

/ / calculation 
f = l/(l + °/ 0 e~(Eg/(2*kT) ) ) ; 

disp(f , ’probability of an e being thermally promoted 
to conduction band is ’); 


Scilab code Exa 18.3 calculation of fraction of e in CB 
clc ; clear ; 
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//Example 18.3 

//calculation of fraction of e in CB 


//given values 
T=300;//temp in K 

kT = . 026 ; //temperture equivalent at room temp in eV 
Egl = . 72;// forbidden gap of germanium in eV 
Eg2 = l . 1;// forbidden gap of silicon in eV 
Eg3 = 5.6;// forbidden gap of diamond in eV 


/ / calculation 

f1 = %e ~(-Egl/(2*kT)) ; 

disp(fl, ’fraction of e in conduction band of 
germanium is’); 

f 2 = %e ~(-Eg2/(2* kT)) ; 

disp(f2, ’ fraction of e in conduction band of 
silicon is’); 

f 3 = %e ~(-Eg3/(2* kT)) ; 

disp(f3, ’fraction of e in conduction band of 
diamond is’); 


Scilab code Exa 18.4 calculation of fractionional change in no of e 

clc ; clear ; 

//Example 18.3 

//calculation of fractionional change in no of e 

//given values 
T1 =300 ; //temp in K 
T2=310 ; //temp in K 

Eg = 1.1; //forbidden gap of silicon in eV 
k = 8.6* 10~-5 ; //boltzmann ’ s constant in eV/K 

/ / calculation 

nl=(10~21.7)*(T1~(3/2))* 10 ~(-2500*Eg/Tl); //no of 
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conduction e at T1 

n2=(10“21.7)* (T2 "(3/2))* 10 ~(-2500*Eg/T2); //no of 
conduction e at T2 
x=n2/nl; 

disp (xtract ional change in no of e is’); 


Scilab code Exa 18.5 calculation of resistivity 

clc ; clear ; 

//Example 18.5 

//calculation of resistivity 

//given values 

e = 1.6*10"-19; 

ni =2.5* 10" 19 ; //i nt r i ns i c density of carriers per in' 3 
ue=. 39; //mobility of e 
uh=. 19; //mobility of hole 


/ / calculation 

c = e*ni*(ue + uh) ; // conductivity 
r=l/c; //resistivity 
d i s p (r , ’resistivity in ohm m is’); 


Scilab code Exa 18.6 calculation of conductivity of intrinsic and doped semiconduc 

clc ; clear ; 

//Example 18.6 

//calculation of conductivity of intrinsic and doped 
semiconductors 

//given values 

h = 4.52* 10"24 ; //no of holes per in"3 


57 




7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

1 

2 


e=1.25*10~14;//no of electrons per m~3 
ue=.38;//e mobility 
uh = . 18 ; / / h o 1 e mobility 


q=l.6*10“-19; //charge of 
/ / calculation 
ni=sqrt(h*e) ; //intrinsic 
ci = q*ni*(ue + uh) ; 
disp (ci , ’ conductivity of 
cp=q*h*uh; 

disp (cp , ’ conductivity of 
) is ’) ; 


e in C 

concentration 

semiconductor ( in S/m) is’); 
doped semiconductor (in S/m 


Scilab code Exa 18.7 calculation of hole concentration 

clc ; clear ; 

//Example 18.7 

//calculation of hole concentration 
//given values 

ni =2.4*10“19; //carrier concentration per m' 3 
N = 4* 10~28 ;// concentration of ge atoms per in' 3 

II calculation 

ND=N/10~6; //donor cocntrtn 
n = ND;//no of elect rones 

p = ni~ 2/n; 

disp(p,’concentartion of holes per m' 3 is’); 


Scilab code Exa 18.8 calculation of Hall voltage 


clc ; clear ; 

//Example 18.8 
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//calculation of Hall voltage 


//given values 

ND = 10~21; //donor density per m' 3 
B=. 5; //magnetic field in T 
J = 500 ; // current density in A/m" 2 
w=3* 10~-3 ; //width in m 
e=l.6*10"-19;// charge in C 

/ / calculation 


V=B*J*w/(ND*e);//in volts 
disp(V*10~3, ’Hall voltage in mv is’); 
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Chapter 19 

PN Junction Diode 


Scilab code Exa 19.1 calculation of potential barrier 

clc ; clear ; 

//Example 19.1 

//calculation of potential barrier 

//given values 

e = 1.6*10"-19; 

n = 4.4* 1 0 ~ 28 ; //no of atoms per m ' 3 

kT= . 026*e ; //temp eqvlnt at room temp 

ni=2.4*10"19; //no of intrinsic carriers per m' 3 

NA=n/10~6;//no of acceptors 

ND=n/10~6; //no of donors 

/ / calculation 

V=(kT/e)*log(NA*ND/ni~2) ; 
d i s p (V , ’ potential barrier in volts is’); 


Scilab code Exa 19.2 calculation of current 
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clc ; clear ; 

//Example 19.2 
//calculation of current 

//given values 

e = 1.6*10"-19; 

kT= . 026*e ; //temp eqvlnt at room temp 

Io = 2*10~-7;// current flowing at room temp in A 

V=. 1;// forward bias voltage in volts 

/ / calculation 

I = Io*(°/„e~(e*V/kT)-l) ; //in Ampere 
disp(I*10~6, ’current flowing when forward bias 
applied (in microampere ) i s ’) ; 
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Chapter 21 

Magnetic Materials 


Scilab code Exa 21.1 calculation of magnetizing force and relative permeability 

clc ; clear ; 

//Example 21.1 

//calculation of magnetizing force and relative 
permeability 

//given values 

M = 2300; //magnetization in A/m 
B= . 00314 ;// flux density in Wb/m'2 
u=12.57*10~-7; //permeability in H/m 

/ / calculation 

H=(B/u)-M; 

di sp (H , ’ m a g net iz ing force (in A/m) is ’); 

Ur=B/(u*H); 

disp(Ur, ’relative permeability is ’) 


Scilab code Exa 21.2 calculation of magnetization and magnetic flux density 
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clc ; clear ; 

//Example 21.2 

//calculation of magnetization and magnetic flux 
density 

//given values 

H=10~5; //external field in A/m 
X=5*10~-5; / / susceptibility 
u=12.57*10~-7; //permeability in H/m 

/ / calculation 

M=X*H; 

disp(M,’ magnetization (in A/m) is ’ ); 

B=u*(M+H); 

disp(B, ’magnetic flux density (in wb/nT2) is ’) 


Scilab code Exa 21.3 calculation of relative permeability 

clc ; clear ; 

// Example 21.3 

/ / calculation of relative permeability 
//given values 

X = 3.7* 10 ~-3 ; //s u s c e p t i b i 1 i t y at 300k 
T=300;//temp in K 
T1 =200 ; //temp in K 
T2 =500 ; //temp in K 

/ / calculation 

C=X*T; //curie constant 

XT1=C/T1; 

disp( XT1, ’relative permeability at T1 is ’); 
XT2=C/T2; 

disp( XT2, ’relative permeability at T2 is’) 
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Chapter 22 

Superconductivity 


Scilab code Exa 22.1 calculation of magnetic field 

clc ; clear ; 

//Example 2 2.1 

//calculation of magnetic field 
//given values 

Tc=7.2; / / transition temp i n K 
T = 5 ; //temp in K 

He =3.3* 10~4 ; //magnetic field at T in A/m 


/ / calculation 

HcO=Hc/(1-(T~2/Tc~2)); 

disp(HcO, ’max value of H at OK (in A/m) is ’); 


Scilab code Exa 22.2 calculation of transition temperature 
clc ; clear ; 
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// Example 2 2.2 

//calculation of transition temperature 
//given values 
T = 8 ; //temp in K 

He = 1 * 10 ~ 5 ; //c r i t i c a 1 magnetic field at T in A/m 
Hc0 = 2*10~5; //magnetic field at 0 K in A/m 

/ / calculation 

Tc=T/ (sqrt (1-Hc/HcO)); 

disp(Tc, ’transition temp in K is’); 


Scilab code Exa 22.3 calculation of temp at which there is max critical field 

clc ; clear ; 

//Example 2 2.3 

//calculation of temp at which there is max 
critical field 

//given values 

Tc=7.26;// critical temp in K 

He =8* 10 ~ 5 ; //max critical magnetic field at T in A/m 
H=4*10~4;// subjected magnetic field at in A/m 

/ / calculation 

T=Tc*(sqrt (1-H/Hc)); 

disp(T, ’max temp for superconductivity in K is’); 
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Chapter 23 
Dielectrics 


Scilab code Exa 23.1 calculation of relative permittivity 

clc ; clear ; 

//Example 2 3.1 

//calculation of relative permittivity 
//given values 

E = 1000; //electric field in V/m 
P=4. 3 * 10 ~ -8 ; //p o 1 arizati o n in C/m'2 
e =8.85*10~-12; //permittivity in F/m 


/ / calculation 

er=l+(P/(e*E)); 

disp(er, ’relative permittivity of NaCl is ’) ; 


Scilab code Exa 23.2 calculation of electronic polarizability 
clc ; clear ; 
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// Example 2 3.2 

/ / calculation of electronic polarizability 
//given values 

e =8.85 *10~-12; //permittivity in F/m 

er = l . 0024; //relative permittivity at NTP 

N = 2.7*10“25; //atoms per m' 3 


/ / calculation 

alpha=e*(er-l)/N; 

disp (alpha electronic polarizability (in F/nP2)is 

); 


Scilab code Exa 23.3 calculation of electronic polarizability and relative permitt 

clc ; clear ; 

//Example 2 3.3 

//calculation of electronic polarizability and 
relative permittivity 

//given values 

e =8.85*10~-12; //permittivity in F/m 
N = 9.8*10~26 ; //atoms per m' 3 
r=.53*10~-10; //radius in m 


/ / calculation 

alpha = 4*°/ 0 pi*e*r~3; 

disp (alpha electronic polarizability (in F/nT2)is 

); 

er = 1+ (4* °/ 0 pi *N*r ~3) ; 

disp(er, ’relative permittivity is ’) 
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Scilab code Exa 23.4 calculation of electronic polarizability and relative perraitt 


clc ; clear ; 

//Example 2 3.4 

//calculation of electronic polarizability and 
relative permittivity 

//given values 

w = 32;// atomic weight of sulphur 
d = 2.08*10~3 ;// density in kg/m'3 
NA=6.02*10~26; //avogadros number 

alpha=3.28*10~-40; //electronic polarizability in F.m 
' 2 

e=8.854*10~-12; //permittiviy 
/ / calculation 

n=NA*d/w; 
k=n*alpha/(3*e); 
er = (l + 2*k)/(l-k) ; 

disp(er, ’relative permittivity is ’) 


Scilab code Exa 23.5 calculation of ionic polarizability 

clc ; clear ; 

//Example 2 3.5 

//calculation of ionic polarizability 

//given values 

n = 1.5; / /refractive index 

er = 6.75; / / relative permittivity 
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/ / calculation 

Pi=(er-n~2)*100/(er-l); 

disp(Pi, ’percentage ionic polarizability (in %)) is’ 

) 


Scilab code Exa 23.6 calculation of frequency and phase difference 

clc ; clear ; 

//Example 2 3.6 

//calculation of frequency and phase difference 

//given values 

t=18*10~-6; //relaxation time in s 

/ / calculation 

f = l/(2*°/ 0 pi*t) ; 

disp(f , ’frequency at which real and imaginary part 
of complx dielectric constant are equal is’); 

alpha = at an (1) * 180/°/„pi ; // phase difference between 
current and voltage ( 1 because real and imaginry 
parts are equal of the dielectric constant) 

di sp (alpha ,’ phase diffeerence (in degree) is’); 


Scilab code Exa 23.7 calculation of frequency 

clc ; clear ; 

//Example 2 3.7 
//calculation of frequency 

//given values 

t=5.5*10~-3; //thickness of plate in m 
Y = 8* 10~ 10 ; //Young ’ s modulus in N/uT2 
d = 2.65*10~3 ;// density in kg/m' 3 
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12 // calculation 

13 f =sqrt (Y/d)/(2*t) ; //in Hz 

14 disp(f/10~3, ’frequency of fundamental note (in KHz) 

is ’) ; 
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Chapter 24 
Fibre optics 


Scilab code Exa 24.1 Fiber optics numerical aperture calculation 

clc ; clear ; 

//Example 24.1 
// Fiber optics 

//given values 

n = 1.5; / /refractive index 

x=.0005; // fractional index difference 

/ / calculation 

u=n*(1-x); 

disp (u cladding index is’); 
alpha=asin(u/n) *180/°/ 0 pi ; 

disp(alpha,’ critical internal reflection angle (in 
degree) is ’); 

theta=asin(sqrt (n~2-u~2))*180/%pi; 

disp (theta, ’ critica 1 acceptance angle(in degree) is’ 

) ; 

N = n* sqrt (2*x) ; 

disp (N ,’numerical aperture is’); 
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Scilab code Exa 24.2 calculation of acceptance angle 

clc ; clear ; 

//Example 24.2 

//calculation of acceptance angle 
//given values 

n = 1.59; //cladding refractive index 
u = l . 33 ;// refract ive index of water 
N = . 20; //numerical aperture offibre 
/ / calculation 

x=sqrt(N~2+n~2); //index of fibre 

Nl = sqrt(x~2-n~2)/u; //numerical aperture when fibre 
is in water 

alpha=asin(Nl) *180/%pi; 

disp (alphaacceptance angle in degree is’); 


Scilab code Exa 24.3 calculation of normalised frequency 

clc ; clear ; 

//Example 24.3 

//calculation of normalised frequency 
//given values 

n=1.45; //core refractive index 

d = . 6 ; // core diametre in m 

N=. 16; //numerical aperture of fibre 

1= . 9*10“-6 ; //wavelength of light 

/ / calculation 

u=sqrt(n~2+N~2);// index of glass fibre 
V = °/ 0 pi *d* sqrt (u~2-n“2)/l; 
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d i s p (V , ’ normalised frequency is’); 


Scilab code Exa 24.4 calculation of normalised frequency and no of modes 

clc ; clear ; 

//Example 24.4 

//calculation of normalised frequency and no of 
modes 

//given values 

n = 1.5 2; //core refractive index 
d=29*10“-6;//core diametre in m 
1 = 1.3*10"-6 ;// wavelength of light 
x=.0007; // fractional refractive index 

/ / calculation 

u=n*(1-x);// index of glass fibre 
V = y„pi *d* sqrt (n“2-u“2)/l; 
d i s p (V , ’normalised frequency is’); 

N=V“2/2; 

disp(N,’no of modes is’); 


Scilab code Exa 24.5 calculation of numerical aperture and maximum acceptance angl 

clc ; clear ; 

//Example 24.5 

//calculation of numerical aperture and maximum 
acceptance angle 

//given values 

n = 1.480; //core refractive index 
u=l.47 ;// index of glass 
1=850*10“-9 ;// wavelength of light 
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V = 2.405 ; //V—number 
/ / calculation 

r = V*l/sqrt (n~2-u~2)/°/„pi/2;//in m 
disp(r*10~6,’core radius in micrometre is’); 
N=sqrt (n~2-u~2); 

disp (N ,’ numerical aperture is’); 
alpha=asin(N) *180/°/„pi; 

di sp (alpha , ’ max acceptance angle is’); 


Scilab code Exa 24.6 calculation of power level 

clc ; clear ; 

//Example 24.6 

//calculation of power level 
//given values 

a = 3.5; //attenuation in dB/krn 

Pi =.5*10~-3; //initial power level in W 

1=4; //length of cable in km 

II calculation 

Po=Pi*10~6/(10~(a*l/10)); 

disp (Po ,’power level after km( in microwatt) is’); 


Scilab code Exa 24.7 calculation of power loss 

clc ; clear ; 

//Example 24.7 
//calculation of power loss 

//given values 

Pi=l*10~-3;// initial power level in W 
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7 1 =. 5; / /length of cable in km 

8 Po=.85*Pi 

9 

10 // calculation 

11 a=(10/l)*logl0(Pi/Po) ; 

12 disp(a,’loss in dB/km is’); 
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Chapter 25 

Digital electronics 


Scilab code Exa 25.1 sura of two binary numbers 

clc ; clear ; 

//Example 2 5.1 

//calculation of sum of two binary numbers 
//given values 

X=’0011 ’;// first binary number 
Y= ’ 0101’; //second binary number 

/ / calculation 

x = bin2dec (X) ; / / decimal equivalent 
y = bin2dec (Y) ; / / decimal equivalent 
z=x+y; 

Z = dec2bin(z) ; 

disp(Z, ’Sum of the given binary numbers is ’) 


check Appendix AP 1 for dependency: 
bin21dec.sci 

check Appendix AP 2 for dependency: 

dec21bin.sci 
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Scilab code Exa 25.2 sura of two binary numbers 


clc ; 
clear ; 

/ / example 2 5.2 

//addition of binary numbers 


a=1010.00; 
b = 00 11.11; 
A=bin2Idee(a); 
B = bin21dec(b) ; 
S=A+B; 


//first number 
//second number 

//converting a in to decimal number 
//converting b in to decimal number 
//adding the two decimal numbers 
temp = dec21bin (S) ; //converting the decimal sum back 
to binary 

f ormat ( ’v ’ , 10) ; //changing the default precision 
to 8 

disp (temp , ’sum is’); //displaying the final output 


check Appendix AP 1 for dependency: 
bin21dec.sci 

check Appendix AP 2 for dependency: 

dec21bin.sci 


Scilab code Exa 25.3 sura of two binary numbers 

clc ; 
clear ; 

/ / example 2 5.3 

//addition of two binary numbers 
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a=1011.01; 
b=1101.11; 
A=bin2Idee(a); 
B = bin21dec(b) ; 
S=A+B; 


// first number 
//second number 

//converting a in to decimal number 
//converting b in to decimal number 
//adding the two decimal numbers 
temp = dec21bin (S) ; //converting the decimal sum back 
to binary 

f ormat ( ’v ’ , 10) ; //changing the default precision 
to 8 

di sp (temp , ’ sum is’); //displaying the final 

output 


Scilab code Exa 25.4 difference of two binary numbers 

clc ; clear ; 

//Example 2 5.4 

//calculation of difference of two binary numbers 
//given values 

X=’1011 ’;// first binary number 
Y= ’ 0101 ’;// second binary number 

/ / calculation 

x = bin2dec (X) ; / / decimal equivalent 
y = bin2dec (Y) ; / / decimal equivalent 
z=x-y; 

Z = dec2bin (z) ; 

disp (Zdifference of the given binary numbers ’); 


Scilab code Exa 25.5 difference of two binary numbers 


clc ; clear ; 

// Example 2 5.5 
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//calculation of difference of two binary numbers 
//given values 

X= ’ 1000’; //first binary number 
Y= ’ 0011 ’;// second binary number 

/ / calculation 

x = bin2dec (X) ; / / decimal equivalent 
y = bin2dec (Y) ; / / decimal equivalent 
z=x-y; 

Z = dec2bin (z) ; 

disp (Z, ’ difference of the given binary numbers ’) 


check Appendix AP 1 for dependency: 
bin21dec.sci 

check Appendix AP 2 for dependency: 

dec21bin.sci 


Scilab code Exa 25.6 difference of two binary numbers 

clc ; 
clear ; 

/ / example 2 5.6 
//binary substraction 

f ormat (’ v 8) ;// changing the default precision to 8 
a = 1001.01;// first number 
b = 00 11 . 10 ;// second number 

A = bin21dec (a) ;// converting a in to decimal number 
B=bin21dec(b);// converting b in to decimal number 
S = A-B; //multiply the two decimal numbers 

temp = dec21bin (S) ;// converting the decimal product 
back to binary 
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disp (temp difference is ’);// displaying the final 
output 


Scilab code Exa 25.7 product of two binary numbers 

clc ; clear ; 

//Example 2 5.7 

//calculation of product of two binary numbers 
//given values 

X=’10101 ’;/ / first binary number with last two digits 
in fractional part 

Y= ’ 101 ’;// second binary number with last two digits 
in fractional part 

/ / calculation 

x = bin2dec (X) ; / / decimal equivalent 
y = bin2dec (Y) ; / / decimal equivalent 
z=x*y; 

Z = dec2bin(z) ; 

disp(Z, ’product of the given binary numbers is ’) 


check Appendix AP 1 for dependency: 
bin21dec.sci 

check Appendix AP 2 for dependency: 

dec21bin.sci 


Scilab code Exa 25.8 binary multiplication 

clc ; 
clear ; 
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/ / example 2 5.8 
//binary multiplication 

f ormat (’v8) ;// changing the default precision to 8 
a = 10101.01;// first number 
b = l 10.10 ;// second number 

A = bin21dec (a) ;// converting a in to decimal number 
B = bin21dec (b) ;// converting b in to decimal number 
S = A*B; //multiply the two decimal numbers 

temp = dec21bin (S) ;// converting the decimal product 
back to binary 

disp (temp product is ’);// displaying the final 
output 


Scilab code Exa 25.9 binary division 

clc ; clear ; 

//Example 2 5.9 

//calculation of quotient of two binary numbers 

//given values 
X= ’1101001 ’ ; // divident 
Y= ’ 101 ’ ; //divisor 

/ / calculation 

x = bin2dec (X) ; / / decimal equivalent 
y = bin2dec (Y) ; / / decimal equivalent 
z=x/y; 

Z = dec2bin(z) ; 

disp(Z, ’quotient of the given binary numbers with 
last two digits in fractional part is ’) 


check Appendix AP 1 for dependency: 
bin21dec.sci 
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check Appendix AP 2 for dependency: 

dec21bin.sci 


Scilab code Exa 25.10 binary division 

clc ; 
clear ; 

//example 2 5.10 
//binary division 

f ormat (’ v8) ;// changing the default precision to 8 
a = 11001 ; // first number 
b = 100 ;// second number 

A = bin21dec (a) ;// converting a in to decimal number 
B = bin21dec (b) ;// converting b in to decimal number 
S=A/B; //multiply the two decimal numbers 

temp = dec21bin (S) ;// converting the decimal product 
back to binary 

d i s p (t emp , ’ quotient is’); //displaying the final 
output 


Scilab code Exa 25.11 octal addition 

clc ; clear ; 

//Example 2 5.11 

//calculation of sum of two octal numbers 

//given values 
X= ’256 ’ ; //divident 
Y= ’437 ’ ; //divisor 
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/ / calculation 

x = oct2dec (X) ; / / decimal equivalent 
y=oct2dec (Y) ; / / decimal equivalent 
z=x+y; 

Z = dec2oct (z) ; //binary equivalent 

disp(Z, ’sum of the given octal numbers is ’) 


Scilab code Exa 25.12 octal multiplication 

clc ; clear ; 

//Example 2 5.12 

//calculation of product of two octal numbers 


//given values 
X=’15’; //divident 
Y= ’ 24 ’ ; //divisor 

/ / calculation 

x = oct2dec (X) ; / / decimal 
y=oct2dec (Y) ; / / decimal 
z=x*y; 

Z=dec2oct(z); //binary 
disp (Zproduct of the 


equivalent 

equivalent 

quivalent 

given octal numbers is ’) 


Scilab code Exa 25.13 hexadecimal addition 

clc ; clear ; 

//Example 2 5.13 

//calculation of sum of hexadecimal numbers 

//given values 

XI = ’C’ ; 

X2= ’A’ ; 
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X3=’E’; 
Y1=’3 ’ ; 
Y2 = ’2 ’ ; 
Y3=’D’ ; 


/ / calculation 

xl=hex2dec (XI) 
x2=hex2dec (X2) 
x3=hex2dec (X3) 
yl=hex2dec (Y1) 
y2=hex2dec (Y2) 
y3=hex2dec (Y3) 
zl=xl+yl 
z2=x2+y2 
z3=x3+y3 
Zl=dec2hex (zl) 
Z2=dec2hex (z2) 
Z3=dec2hex (z3) 
disp (Zl , ’sum of 
is’); 

disp (Z2 , ’ sum of 
numbers is’); 
di sp (Z3 , ’ sum of 
numbers is’); 


// decimal 
// decimal 
// decimal 
// decimal 
// decimal 
// decimal 


equivalent 

equivalent 

equivalent 

equivalent 

equivalent 

equivalent 


//binary equivalent of sum 
//binary equivalent of sum 
//binary equivalent of sum 
the first set of hexadecimal numbers 

the second set of hexadecimal 

the t h i r d m set of hexadecimal 


Scilab code Exa 25.14 binary to decimal conversion 

clc ; clear ; 

//Example 2 5.13 

//conversion of binary to decimal 

//given values 

X=10.101;// binary number 

II calculation 
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Z=(l*2~l) + (0*2~0)+(1*2~-1)+(0*2~-2)+(l*2~-3) ; 
d i s p (Z , ’ decimal equivalent of the given binary 
number is’) 


Scilab code Exa 25.15 decimal to binary conversion 

clc ; clear ; 

//Example 2 5.15 

//conversion of decimal to binary 

//given values 
X=43; //decimal number 

/ / calculation 

Z = dec2bin(X) ; 

d i s p (Z , ’ binary equivalent of the given decimal 
number is’); 


Scilab code Exa 25.16 decimal to binary conversion 

clc;//clears the command window 
clear ;// clears all the variables 
//example 2 5.16 

//decimal to binary conversion 

format ( ’v ’ , 18); //changing the default prcision to 20 
significant digits 


i = l;x = l; // flag bits 

dec =43.3125 ;// given decimal number which should be 
expressed in binary 
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temp2 = f loor (dec) ;// separating integer part from the 
given number 

temp4 = modulo (dec , 1) ;// separating decimal part from 
the given number 

while (terap2>0) //storing each integer digit in vector 
for convenience 

p(i) = ( modulo ( floor (temp2) ,2)) 
temp2=f loor (temp2/2); 

1 = 1 + 1 ; 

end 

temp2=0;// clearing temporary variable ’temp2’ 

for j=1: length (p) 

//multipliying bits of integer part with their 
position values and adding 

temp2=temp2+(p(j)*10~(j-l)); 

end 

while(temp4~ = 0) //storing each decimal digit in 
vector for convenience 

temp4=temp4*2; 
d(x)=floor(temp4) ; 
x=x+l; 

temp4 = modulo (temp4 , 1) ; 

end 

temp5=0;// clearing temporary variable ’temp5’ 

for j=1: length (d) 

//multipliying bits of decimal part with their 
position values and adding 

temp5 = temp5 + (10~(-1 * j)*d(j)) 

end 

temp3=temp2+temp5; 

//finally adding both the integer and decimal parts 


87 


43 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 


to get total output. 

disp(temp3, ’the equivalent binary number is ’); 


Scilab code Exa 25.17 decimal to octal conversion 

clc;// clears the command window 
clear ;// clears all the variables 
/ / example 2 5.17 
//decimal to octa conversion 

f ormat (’v8) ;// making the default precision to 8 
significant digits 

i=l;w=l; 

dec =375.23 ;// given decimal number which should be 
expressed in base 8 

temp = modulo (dec , 1) ;// separating decimal part from 
the given number 

temp2 = floor (dec) ;// separating integer part from the 
given number 


while (temp2>0) //storing each integer digit in vector 
for convenience 

p(i) = ( modulo ( floor (temp2) ,8)) 
temp2 = floor (temp2/8) ; 
i=i+l; 

end 

temp2=0 ;// clear ing temporary variable ’temp2’ 

for j=1: length (p) 

//multipliying bits of integer part with their 
position values and adding 

temp2=temp2+(p(j)*10~(j-l)); 


end 
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while (temp ~ = 0) //storing each decimal digit in 
vector for convenience 

terap = temp *8; 
q(w)=floor(temp) ; 
w=w+l; 

terap = modulo (temp,1) ; 

end 

templ=0; //flag bit 

for k=l: length (q) 

//multipliying bits of decimal part with their 
position values and adding 

templ=templ+(10~(-l*k)*q(k)); 

end 

temp3=temp2+templ; 

disp (temp3 octal number is’); 


Scilab code Exa 25.18 octal to binary conversion 

clc ; clear ; 

//Example 2 5.18 
// oca11 to binary conversion 

//given values 
X=’257’;//octal number 

/ / calculation 

x = oct2dec (X) ; / / decimal equivalent 
Z = dec2bin(x) ; 

disp(Z, ’binary number is ’) 


Scilab code Exa 25.19 octal to binary conversion 
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clc;//clears the command window 
clear ;// clears all the variables 
/ / example 2 5.19 
//octal to binary conversion 

format ( ’ v ’ , 8); //setting the default precision to 8 

i=l;w=l; 

bin = 34.56; //Given octal number which we need to be 
convert into binary 

templ = f loor (bin) ; // separating integer part from the 
given number 

tempO=modulo (bin,1);// separating decimal part from 
the given number 

temp2 = temp0 *10~2 ;// converting decimal value to 
interger for convenience 
while (tempi>0) //storing each integer digit in 

vector for convenience 

p(i) =modulo (tempi ,10) ; 
templ=round (tempi/10); 
i=i+l; 

end 

while (temp2 >0) //storing each decimal digit in 

vector for convenience 

q(w)=modulo(temp2 ,10) ; 
temp2 = floor (temp2/10) ; 
w=w+1; 

end 

templ=0; //clearing temporary variable ’tempi 

for i=1:2 

//multipliying bits of decimal part with their 
position values and adding 

templ = templ + (p(i)*8~(i-l)) ; 

end 
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temp2=0 ;// clearing temporary variable ’ temp2 ’ 

for z = 1:2 

//multipliying bits of decimal part with their 
position values and adding 

temp2 = temp2 + (q(z)*8~(-l*(3-z))) ; 


end 

temp=templ+temp2; 

//adding both integer and decimal parts to get total 
deciaml value . 

dec=temp; 

t emp2 = f loor (dec ) ; //separating integer part from the 
given number 

temp3 = modulo (dec , 1) ;// separating decimal part from 
the given number 

format ( ’v 18) ;// setting the default precision to 8 
i = l; x = l ; //flag bits 

while (terap2>0) //storing each integer digit in vector 
for convenience 

p(i) = ( modulo ( floor (temp2) ,2)) 
temp2=f loor (temp2/2); 

1 = 1 + 1 ; 

end 

temp2=0; //clears temporary variable ’temp2 ’ 

for j=1: length (p) 

//multipliying bits of integer part with their 
position values and adding 

temp2=temp2+(p(j)*10~(j-l)); 

end 

temp4 = modulo (temp3 , 1) ; 
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wh ile(te rap4~ = 0) //storing each decimal digit in 
vector for convenience 

temp4=temp4*2; 
d(x)=floor(temp4) ; 
x=x+l; 

temp4 = modulo (temp4 , 1) ; 

end 

temp5=0; //clears temporary variable ’temp2 ’ 

for j=1: length (d) 

//multipliying bits of decimal part with their 
position values and adding 

temp5=temp5+(10~(-1 * j)*d(j)) 

end 

temp=temp2+temp5; 

//finally adding both the integer and decimal parts 
to get total output, 
disp (temp binary number is’); 


Scilab code Exa 25.20 binary to octal conversion 

clc;// clears the command window 

clear ;// clears all the variables 

i=l;w=l; 

bin = 1011.01101 ;// Given binary number which we need 
to be convert into octal 

//conversion to decimal first 

templ = floor (bin) ;// separating integer part from the 
given number 

temp2 = modulo (bin , 1) ;// separating decimal part from 
the given number 
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temp2 = temp2 * 10~5 ;// converting decimal value to 
integer for convenience 

while (tempi>0) //storing each integer digit in vector 
for convenience 

p(i) =modulo (tempi ,10) ; 
templ = floor (tempi/10) ; 
i=i+l; 

end 

while (temp2 >0) //storing each decimal digit in vector 
for convenience 

q(w)=modulo(temp2 ,2) ; 
temp2=(temp2/10); 
temp2=f loor (temp2); 
w=w+l; 

end 

tempi =0 ;// flag bit 

for i = l: length (p) //checking whether it is a binary 
number or not 
if (p (i) >1) then 

disp(’not a binary number’); 
abort ; 

end 

end 

for i=l: length (p) 

//mnltipliying bits of integer part with their 
position values and adding 

templ = templ + (p(i)*2~(i-l)) ; 

end 

temp2=0; //f1ag bit 

for z=l: length (q) 

//multipliying bits of decimal part with their 
position values and adding 

temp2 = temp2 + (q(z)*2~(-l*(6-z))) ; 

end 

dec=templ+temp2; 

//finally adding both the integer and decimal parts 
to get decimal equivalent 
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//conversion from decimal to octal 

f ormat (’ v 8) ;// making the default precision to 8 
significant digits 

i=l;w=l; 

temp = modulo (dec , 1) ;// separating decimal part from 
the given number 

temp2 = f loor (dec) ;// separating integer part from the 
given number 


while (terap2>0) //storing each integer digit in vector 
for convenience 

r(i) = ( modulo ( floor (temp2) ,8)) 
temp2 = floor (temp2/8) ; 
i=i+l; 

end 

temp2=0;// clearing temporary variable ’temp2’ 

for j=1: length (r) 

//multipliying bits of integer part with their 
position values and adding 

temp2=temp2+(r(j)*10~(j-l)); 

end 

while (temp ~ = 0) //storing each decimal digit in 
vector for convenience 

temp = temp *8; 
s(w)=floor(temp) ; 
w=w+1; 

temp=modulo (temp,1); 

end 

templ=0; //flag bit 

for k=l: length (s) 

//multipliying bits of decimal part with their 
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position values and adding 

templ = templ + (10~(-l*k)*s(k)) ; 

end 

temp3 = terap2 + templ ; 

disp (temp3 octal number is’); 


Scilab code Exa 25.21 hexa to decimal conversion 

clc ; clear ; 

//Example 2 5.21 

//hexadecimal to decimal conversion 
//given values 

X= ’AC5 ’ ; // hexadecimal number 
/ / calculation 

x = hex2dec (X) ; / / decimal equivalent 
disp(x, ’decimal number is ’) 


Scilab code Exa 25.22 decimal to hexadecimal conversion 

clc;//clears the command window 
clear ;// clears all the variables 
//example 2 5.22 

//decimal to hexadecimal conversion 
f ormat (’v4) ;// making the default precision to 8 
significant digits 
dec =379 .54 ;// given decimal 
w=l;i=1; 

templ = floor (dec) ;// separating integer part from the 
given number 
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temp2 = modulo (dec , 1) ;// separating decimal part from 
the given number 

x=dec2hex (tempi);// hexadecimal equivalent of integer 
part 

s =0; 

while(terap2~ = 0) //storing each decimal digit in 
vector for convenience 

temp2=temp2*16; 
q(w)=floor(temp2) ; 
s = s +1; //counter of a 
a(w)=dec2hex(q(w)) ; 
w=w+l; 

temp2 = modulo (temp2 ,1) ; 

end 

f=a(1); 

for i=2:s 

f = f + a(i) ; 

end 

b=’.’;// for concatenating to get the decimal part of 
hexadecimal 

hex = x + b + f ;// concatenating integer and decimal part 
disp (hex hexadecimal equivalent is’); 


Scilab code Exa 25.23 hexa to binary conversion 

clc ; clear ; 

//Example 2 5.23 

//hexadecimal to binary conversion 
//given values 

X=’7AB’ ; //hexadecimal number 
/ / calculation 

x = hex2dec (X) ; / / decimal equivalent 
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z = dec2bin(x) ; 

disp(z, ’binary number is ’); 


Scilab code Exa 25.24 binary to hexa conversion 

clc ; clear ; 

//Example 2 5.24 

//binary to hexadecimal conversion 

//given values 

X= ’ 1011101 ’ ; // binary number 

/ / calculation 

x = bin2dec (X) ; / / decimal equivalent 
z = dec2hex (x) ; 

disp(z, ’hexadecimal number is ’); 


Scilab code Exa 25.25 Substraction by ones complement method 


clc ; 
clear ; 

//Example 2 5.25 

//substraction by one’s complement method 
//aaa=input (” Enter the first no (in decimal) 
//bb=input(” Enter the number from which first no 
has to be substracted:”); 

bb=14; 

aaa = -7;// substraction is addition of negative 
number 

if aaa<0 then 
aa=-l*aaa; 
else aa=aaa; 
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//converting from decimal to 


end 
a = 0; 
b = 0; 
q = 0 ; 

for i=l:5 

binary 

x = modulo (aa , 2) ; 
a= a + (10“q)*x; 
aa=aa/2; 
aa=f loor (aa); 
q=q+l; 

end 
q = 0; 

for i = l:5 //converting from decimal to binary 

y = modulo (bb , 2) ; 
b= b + (10“q)*y; 
bb=bb/2; 
bb=f loor (bb); 
q=q+l; 

end 

for i=l:5 

al(i)=modulo(a,10) ; 
a=a/10; 
a=round (a); 


end 

for i=l:5 

bl (i)=modulo (b,10); 
b=b/10; 
b = round (b) ; 

end ; 

if aaa<0 then// making one’s complement if number 
less than zero 

for i=l:5 

al(i)=bitcmp(al(i) ,1) ; 

end 

car(1)=0; 


i s 
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for i=l:5 

cl(i)=al(i)+bl(i) + car(i) ; 
if cl(i)== 2 then 
car(i+l)= 1; 
cl(i)=0; 

elseif cl(i)==3 then 
car(i + 1)= 1; 
cl(i)=1; 

else 

car(i+1)=0; 

end ; 

end ; 

car2(l)=car(6) ; 
re =0; 

format ( ’v ’ , 18) ; 
for i=l:5 

re=re+(cl(i)*(10“(i-l))) //result of one’s 
complement addition 

end ; 


for i=l:5 

s(i)=modulo(re ,10) ; 
re=re/10; 
re=round (re); 

end ; 

for i=l:5 

rel(i)=s(i) + car2(i); //addition of carry after one’s 
complement addition 

if rel(i)== 2 then 
car2(i + 1)= 1; 
re 1(i)=0; 

elseif rel(i)==3 then 
car2(i + 1)= 1; 
rel(i)=1; 
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else 


car2(i+1)=0; 

end ; 

end ; 
re2 =0; 

format ( ’ v ’ ,18) ; 
for i=l:5 

re2=re2+(rel(i)*(10~(i-l))) 

end ; 

disp(re, ’difference is ’) 


Scilab code Exa 25.26 Substraction by ones complement method 


clc ; 
clear ; 

/ / example 2 5.26 

//substraction by one’s complement method 
//a=input(” Enter the first no (binary) :”) ; 
//b=input(” Enter the number from which first no has 
to be substracted ; 

a=10001; 
b=10011; 
q = 0 ; 

for i=l:5 

al(i)=modulo(a,10) ; 
a=a/10; 
a=round (a); 


end 

for i=l:5 

bl (i)=modulo (b,10); 
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b=b/10; 
b = round (b) ; 

end ; 

for i = 1 :5 //making one’s complement of number to be 
substracted 

al(i)=bitcmp(al(i) ,1) ; 

end 

car(1)=0; 
for i=1:5 

cl(i)=al(i)+bl(i)+car(i); 
if c1(i)== 2 then 
car(i+l)= 1; 
cl(i)=0; 

elseif cl(i)==3 then 
car(i + 1)= 1; 
cl(i)=1 ; 

else 

car(i+1)=0; 

end ; 

end ; 

car2(l)=car(6); 
re =0; 

format ( ’v ’ , 18) ; 
for i=1:5 

re=re+(cl(i)*(10~(i-1))) //result of one’s 
complement addition 

end ; 


for i=l:5 

s(i)=modulo(re ,10) ; 
re=re/10; 
re=round (re); 

end ; 

if car2(l)==l then// checking carry 
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for i=l:5 

rel(i)=s(i) + car2(i);// addition of carry after 
complement addition 

if rel(i)== 2 then 
car2(i + 1)= 1; 
re 1(i)=0; 

elseif rel(i)==3 then 
car2(i+l)= 1; 
rel(i)=1; 

else 

car2(i+1)=0; 

end ; 

end ; 


re2 =0; 

format ( ’v ’ , 18) ; 
for i=l:5 

re2=re2+(rel(i)*(10~(i-l))) 

end ; 

disp(re2, ’difference is ’) ; 

else 

for i=l:5 

rel(i)=bitcmp(s(i),1); 

end 

re2 =0; 
for i = l : 5 

re2=re2+(rel(i)*(10~(i-l))) 

end ; 

re2=-l*re2; 

disp(re2, ’difference is ’) ; 

end ; 


Scilab code Exa 25.27 Substraction by ones complement method 
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clc ; 
clear ; 

//example 2 5.27 

//substraction by one’s complement method 
//a=input (” Enter the first no (binary) ; 
//b=input(” Enter the number from which first no has 
to be substracted ; 

a=10011; 
b=10001; 
q = 0 ; 

for i=1:5 

al(i)=modulo(a,10) ; 
a=a/10; 
a=round (a); 

end 

for i=1:5 

bl (i)=modulo (b , 10) ; 
b=b/10; 
b = round (b) ; 

end ; 

for i = 1 : 5 //making one’s complement of number to be 
substracted 

al(i)=bitcmp(al(i) ,1) ; 

end 

car(1)=0; 
for i=l:5 

cl(i)=al(i)+bl(i)+car(i); 
if c1(i)== 2 then 
car(i+l)= 1; 
cl(i)=0; 

elseif cl(i)==3 then 
car(i + 1)= 1; 
cl(i)=1 ; 

else 
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car(i+l)=0; 


end ; 

end ; 

car2(l)=car(6); 
re =0; 

format ( ’v ’ , 18) ; 
for i=l:5 

re=re+(cl(i)*(10“(i-l))) //result of one’s 
complement addition 

end ; 

for i=1:5 

s(i)=modulo(re ,10) ; 
re=re/10; 
re=round (re); 

end ; 

if car2(l)==l then// checking carry 
for i = 1 : 5 

rel(i)=s(i) + car2(i);// addition of carry after one’s 
complement addition 

if rel(i)== 2 then 
car2(i + 1)= 1; 
re 1(i)=0; 

elseif rel(i)==3 then 
car2(i + 1)= 1; 
rel(i)=1; 

else 

car2(i+1)=0; 

end ; 

end ; 


re2 =0; 

format ( ’v ’ , 18) ; 
for i=l:5 

re2=re2+(rel(i)*(10~(i-l))) 

end ; 

re2= -1* re2; 
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disp(re2, ’difference is ’) ; 


else 

for i=l:5 

rel(i)=bitcmp(s(i) ,1); 

end 

re2 =0; 
for i = l : 5 

re2=re2+(rel(i)*(10~(i-l))) 

end ; 

re2=-l*re2; 

disp(re2, ’difference is ’) ; 


end ; 


Scilab code Exa 25.28 finding twos complement 

clc ; 
clear ; 

// example25.2 8 

//f ini ding two’s complement 

//a=input(” Enter the number (binary) :”) ; 

a = 1010 ; 
for i=l:4 

al(i)=modulo(a,10) ; 
a=a/10; 
a=round (a); 

end 

for i = 1 : 4 //making one’s complement of number 

al(i)=bitcmp(al(i) ,1); 

end 

for i=l:4 
car(1)=1; 

re ( i) = al (i) + car ( i ) ; //addit io n of one to one’s 
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complement to contain two’s complement 

if re(i)== 2 then 
car(i+l)= 1; 
re(i)=0; 

elseif re(i)==3 then 
car(i + 1)= 1; 
re(i)=1 ; 

else 

car(i+1)=0; 

end ; 

end ; 


re2 =0; 

format ( ’v ’ , 18) ; 
for i=l:4 

re2=re2+(re(i)*(10~(i-l))) 

end ; 

disp (re2, ’two s complement is’); 


Scilab code Exa 25.29 Addition of negative number by twos complement method 


clc ; 
clear ; 

//Example 2 5.29 

//addition of negative number by by two’s complement 
method 

//bb=input(” Enter the first no (in decimal) : ” ) ; 

//aaa=input (” Enter the negative number that has to 
be added”) ; 

bb = 14; 
aaa=-7; 
if aaa<0 then 
aa = -l* aaa; 
else aa=aaa; 
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//converting from decimal to 


end 
a = 0; 
b = 0; 
q = 0 ; 

for i=l:5 

binary 

x = modulo (aa , 2) ; 
a= a + (10“q)*x; 
aa=aa/2; 
aa=f loor (aa); 
q=q+l; 

end 
q = 0; 

for i = l:5 //converting from decimal to binary 

y = modulo (bb , 2) ; 
b= b + (10“q)*y; 
bb=bb/2; 
bb=f loor (bb); 
q=q+l; 

end 

for i=l:5 

al(i)=modulo(a,10) ; 
a=a/10; 
a=round (a); 


end 

for i=l:5 

bl (i)=modulo (b,10); 
b=b/10; 
b = round (b) ; 

end ; 

if aaa<0 then// making one’s complement of the 
negative number 

for i=l:5 

al(i)=bitcmp(al(i) ,1) ; 

end 

car(1)=0; 
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for i=l:5 

cl(i)=al(i)+bl(i) + car(i) ; 
if cl(i)== 2 then 
car(i+l)= 1; 
cl(i)=0; 

elseif cl(i)==3 then 
car(i + 1)= 1; 
cl(i)=1; 

else 

car(i+1)=0; 

end ; 

end ; 
re =0; 

format ( ’v ’ , 18) ; 
for i=l:5 

re=re+(cl(i)*(10“(i-l))) //result of one’s 
complement addition 

end ; 

for i=l:5 

s(i)=modulo(re ,10) ; 
re=re/10; 
re=round (re); 

end ; 

if car(6)==l then// checking carry 
car2(1)=1; 

for i=1:5 

rel(i)=s(i) + car2(i); //addition of carry after one’s 
complement addition 

if rel(i)== 2 then 
car2(i + 1)= 1; 
re 1(i)=0; 

elseif rel(i)==3 then 
car2(i+1)= 1; 
rel(i)=1; 

else 

car2(i+1)=0; 
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end ; 


end ; 


re2 =0; 

format ( ’v ’ , 18) ; 
for i=l:5 

re2=re2+(rel(i)*(10'(i-l))) 

end ; 

disp(re2, ’difference is ’) ; 

else 

for i=l:5 

rel(i)=bitcmp(s(i),1); 

end 

re2 =0; 
for i = l : 5 

re2=re2+(rel(i)*(10~(i-l))) 

end ; 

re2=-l*re2; 

disp(re2, ’difference is ’) ; 

end ; 


Scilab code Exa 25.30 Substraction by twos complement method 


clc ; 
clear ; 

/ / example 2 5.27 

//substarction by one’s complement method 
//a=input(” Enter the first no (binary) :”) ; 
//b=input(” Enter the number from which first no has 
to be substracted :”) ; 

a=10011; 
b=10001; 
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q = 0 ; 


for i=1:5 

al(i)=modulo(a,10) ; 
a=a/10; 
a=round (a); 

end 

for i=l:5 

bl (i)=modulo (b,10); 
b=b/10; 
b = round (b) ; 

end ; 

for i = 1 :5 //making one’s complement of number to be 
substracted 

al(i)=bitcmp(al(i) ,1) ; 

end 

car(1)=0; 
for i=l:5 

cl(i)=al(i)+bl(i)+car(i); 
if c1(i)== 2 then 
car(i+1)= 1; 
cl(i)=0; 

elseif cl(i)==3 then 
car(i + 1)= 1; 
cl(i)=1 ; 

else 

car(i+1)=0; 

end ; 

end ; 
re =0; 

format ( ’v ’ , 18) ; 
for i=l:5 

re=re+(cl(i)*(10“(i-1))) //result of one’s 
complement addition 

end ; 
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for i=1:5 

s(i)=modulo(re ,10) ; 
re=re/10; 
re=round (re); 

end ; 

if car(6)==l then// checking carry 
for i=1:5 

rel(i)=s(i) + car2(i);// addition of carry after 
complement addition 

if rel(i)== 2 then 
car2(i + 1)= 1; 
re 1(i)=0; 

elseif rel(i)==3 then 
car2(i + 1)= 1; 
rel(i)=1; 

else 

car2(i+1)=0; 

end ; 

end ; 


re2 =0; 

format ( ’v ’ , 18) ; 
for i=l:5 

re2=re2+(rel(i)*(10~(i-l))) 

end ; 

re2= -1* re2; 

disp(re2, ’difference is ’) ; 

else 

for i=l:5 

rel(i)=bitcmp(s(i),1); 

end 

re2 =0; 
for i = l : 5 

re2=re2+(rel(i)* (10 ~(i-1))) 

end ; 


one ’ s 
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83 re2 = -l*re2; 

84 disp(re2,’ difference is’) 

85 

86 end ; 
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Appendix 


Scilab code AP 1 Binary to Decimal convertor 

1 //bin21dec is a function whcih converts any binary 

number given to it will output its equivalent 
decimal number 

2 //pass the binary number as an argument to the 

function 

3 // For eg : bin2 1 decimal (10 10) 

4 //Will give an output of 10 

5 

6 function [temp]=bin21dec(bin) 

7 i=l;w=l; 


8 


9 


templ=f loor (bin); 


// separating 


10 


integer part from the given number 

temp2 = modulo (bin , 1) ; 


// separating 


11 


decimal part from the given number 

temp2=temp2*10~3; 


/ / converting 


decimal value to interger for convenience 


12 

13 


while (tempi>0) 


14 

15 


//storing each integer digit in vector for 
convenience 

p(i) =modulo (tempi ,10) ; 
templ = floor (tempi/10) ; 
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i=i+l; 

end 

while (temp2 >0) 

//storing each integer digit in vector for 
convenience 

q(w)=modulo (temp2 , 2) ; 
temp2 = (temp2/10) ; 
t emp2=f loor (temp2); 
w=w+1; 

end 

tempi=0; 

//clearing the temporary variable ’temp2 ’ 

for i=l: length (p) 

//checking whether it is binary or not. 

if (p (i) >1) then 

dispC’not a binary number’); 
abort ; 

end 

end 

for i=l: length (p) 

//mult ipliy ing the bits of integer part with 
their position values and adding 

templ = templ + (p(i)*2~(i-l)) ; 

end 

temp2 =0; 

//clearing the temporary variable ’temp2 ’ 

for z=l:w-1 

//multipliying the bits of decimal part with 
their position values and adding 

temp2 = temp2 + (q(z)*2 ~(-l*(4-z))) ; 

end 


114 


45 temp = templ + temp2 ; 

//finally adding both the integer and decimal 
parts to get total output. 

46 endfunction 


Scilab code AP 2 Decimal to Base 2 Converter 


1 //dec21bin is a function whcih converts any decimal 

number given to it will output its equivalent 
binary number 

2 //pass the decimal number as an argument to the 

function 

3 // For eg:dec21bin(10) 

4 //Will give an output of 1010 

5 


6 function [temp]=dec2lbin(dec) 

7 temp2=f loor (dec); 

integer part from the given 

8 temp4 = modulo (dec , 1) ; 

decimal part from the given 

9 

10 format ( ’v ’ ,18) ; 

the default precision to 18 

11 

12 i = l; p = 0 ; x = l; 

bits 


// separating 
number 

// separating 
number 


// changing 


// flag 


13 

14 while (temp2>0) 

//storing each integer digit in vector for 
convenience 

15 p(i)=( modulo ( floor (temp2),2)) 

16 temp2 = floor (temp2)/2 ; 

17 i = i + l; 

18 end 
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temp2 =0; 

//clearing the temporary variable ’temp2 ’ 

for j=1: length (p) 

//multipliying the bits of integer part with 
their position values and adding 

temp2=temp2+(p(j)*10~(j-l)); 

end 

while (terap4~=0) 

//storing 

each integer digit in vector for convenience 

temp4=temp4*2; 
d(x)=floor(temp4) ; 
x=x+1; 

temp4 = modulo (temp4 , 1) ; 

end 

temp5 =0; 

//clearing the temporary variable ’temp2’ 

for j =1:x-1 // 

multipliying the bits of decimal part with 
their position values and adding 

temp5 = terap5 + (10~( — 1 * j)*d(j)) 

end 

temp=temp2+temp5; 

// 

finally adding both the integer and decimal 
parts to get total output. 

endfunction 
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